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Abstract 

The aim of this paper is to give a formulation of the dynamics of 
nonlinear RLC circuits as a geometric Birkhoffian system and to discuss 
in this context the concepts of regularity, conservativeness, dissipative- 
ness. An RLC circuit, with no assumptions placed on its topology, will 
be described by a family of Birkhoffian systems, parameterized by a fi- 
nite number of real constants which correspond to initial values of certain 
state variables of the circuit. The configuration space and a special Pfaf- 
fian form, called Birkhoffian, are obtained from the constitutive relations 
of the resistors, inductors and capacitors involved and from Kirchhoff 's 
laws. Under certain assumptions on the voltage-current characteristic for 
resistors, it is shown that a Birkhoffian system associated to an RLC cir- 
cuit is dissipative. For RLC networks which contain a number of pure 
capacitor loops or pure resistor loops the Birkhoffian associated is never 
regular. A procedure to reduce the original configuration space to a lower 
dimensional one, thereby regularizing the Birkhoffian, it is also presented. 
In order to illustrate the results, specific examples are discussed in detail. 

Keywords: geometric methods in differential equations, Birkhoffian 
differential systems, Birkhoffian vector fields, dissipative dynamical sys- 
tems, electrical networks 
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1 Introduction 

Lagrangian and Hamiltonian mechanics continue to attract a large amount of 
attention in the literature, because many mechanical and electro-mechanical 
systems may be modelled within these frameworks. During the past ten years, 
a far-reaching generalization of the Hamiltonian framework has been developed 
in a series of papers. This generalization, which is based on the geometric notion 
of generalized Dirac structure, gives rise to implicit Hamiltonian systems (see 
for example papers by Maschke and van der Schaft JU|> van der Schaft |12|L 
In the papers by Yoshimura and Marsden ^2] the concept of Dirac structures 
and variational principle are used to define and develop the basic properties of 
implicit Lagrangian systems. 
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An alternative approach to the study of dynamical systems which appears to 
cover a wide class of systems, among them the nonholonomic systems, the de- 
generate systems and the dissipative ones, is the Birkhoffian formalism, a global 
formalism of the dynamics of implicit systems of second order ordinary differen- 
tial equations on a manifold. The classical book by G. D. Birkhoff [2], contains 
in Chapter I many interesting ideas about classical dynamics from the viewpoint 
of differential geometry. In order to present these ideas in a coordinate free fash- 
ion, one considers the formalism of 2-jets (see for example Kobayashi and Oliva 
0). The space of configurations is a smooth m-dimensional diffcrcntiable con- 
nected manifold and the covariant character of the Birkhoff generalized forces 
is obtained by introducing the notion of elementary work, called Birkhofnan, 
a special Pfaffian form defined on the 2-jets manifold. The dynamical system 
associated to this Pfaffian form is a subset of the 2-jets manifold which defines 
an implicit second order ordinary differential system. The notion of Birkhofnan 
allows to formulate the concepts of reciprocity, regularity, affine structure in the 
accelerations, conservativeness, in an intrinsic way. 

The electrical circuits theory benefits from many tools developed in math- 
ematics. In order to study the dynamics of LC and RLC electrical circuits, 
various Lagrangian and Hamiltonian formulations have been considered in the 
literature (see for example [3], @|, [HI, [TJ]|, |T5], and the references therein). 
To describe LC circuit equations, Hamiltonian formulations have been used more 
often. In the dynamics of a nonlinear LC circuit is shown to be of Hamilto- 
nian nature with respect to a certain Poisson bracket which may be degenerate, 
that is, non-symplectic. The dynamics of "complete" RLC networks was de- 
scribed by Brayton and Moser 0] in terms of a function of inductor currents 
and capacitor voltages, called the mixed potential function. However, for all 
those formulations, a certain topological assumption on the electrical circuit 
appears to be crucial, that is, the circuit is supposed to contain neither loops of 
capacitors nor cutsets of inductors. In [TJJ, an d the Poisson bracket is 
replaced by the more general notion of a Dirac structure on a vector space, lead- 
ing to implicit Hamiltonian systems. In this formalism, it is possible to include 
LC networks which do not obey the topological assumption mentioned before. 
In |15| an example of LC circuit in the context of implicit Lagrangian systems 
is given for a degenerate Lagrangian system with holonomic constraints. 
The potential relevance of the Birkhofnan formalism in the context of electrical 
circuits is discussed by Ionescu and Scheurle [B], where a formulation of general 
nonlinear LC circuits within the framework of Birkhoffian dynamical systems 
on manifolds is presented. In |S] specific examples of electrical networks are 
discussed in this framework. These are networks which contain closed loops 
formed by capacitors, as well as inductor cutsets, and also LC networks which 
contain independent voltage sources as well as independent current sources. 

In the paper at hand we present a formulation of the dynamics of nonlinear 
RLC electrical circuits within the framework of Birkhoffian systems. Based on 
Kirchhoff 's laws and the constitutive relations for the resistors, inductors, capac- 
itors involved, we get for a nonlinear RLC circuit, a whole family of configura- 
tion spaces and a special Pfaffian forms, called Birkhoffians. The configurations 
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spaces are parameterized by a finite number of real constants which correspond 
to initial values of certain state variables of the circuit. In particular, we can 
allow pure capacitor loops as well as pure inductor cutsets. Under certain as- 
sumptions on the voltage-current characteristic for resistors, it is shown that a 
Birkhoffian system associated to an RLC circuit is dissipative. 

The paper is organized as follows. In Section 2, we recall the basics of 
Birkhoffian systems (see presented from the viewpoint of differential ge- 
ometry using the formalism of jets (see [7j). In particular, we introduce the 
notion of dissipative Birkhoffian system, in order to be able to treat the case 
of RLC networks in the next section. In Section 3, our Birkhoffian formulation 
of the dynamic equations of a nonlinear RLC circuit is given. Properties of the 
corresponding Birkhoffian such as its regularity and its dissipativness are also 
discussed in this section. For electrical RLC networks which contain a number 
of only capacitor loops or only resistors loops, we present a systematic procedure 
to reduce the original configuration space to a lower dimensional one, thereby 
regularizing the Birkhoffian. On the reduced configuration space the reduced 
Birkhoffian will still be dissipative. Finally, in section 4 we consider two specific 
examples. These examples are supposed to serve our purpose of demonstrating 
the power of the Birkhoffian approach. 

Acknowledgement. I would like to express my gratitude to Professor Mars- 
den for valuable discussions concerning the topic of this paper. 

2 Birkhoffian systems 

For a smooth m-dimensional differentiable connected manifold M, we consider 
the tangent bundles (TM,tt m ,M) and (TTM,tt tm ,TM). Let q = (q 1 , q 2 ,..., 
q m ) be a local coordinate system on M . This induces natural local coordinate 
systems on TM and TTM, denoted by (q, q), respectively (q, q, dq, dq). The 
2- jets manifold J 2 (M) is a 3m-dimensional submanifold of TTM defined by 

J 2 (M) = {ze TTM | Tir M (z) = tttm(z)} (2.1) 

where T-km ■ TTM — » TM is the tangent map of ttm- We write ttj :— 
ttm|j2(m) = Tttm\ j 2 (M)- (J 2 (M), 7Tj, TM), called the 2-jet bundle (see [7]), is 
an affine bundle modelled on the vertical vector bundle (V(M), ktm\v(m)i TM), 
V(M) = \J veTM V V (M) , where V v (M) = {z e T V TM \ (Ttt m ) v (z) = 0}. In [T], 
this bundle is denoted by T 2 (M) and named second-order tangent bundle. 
In natural local coordinates, the equality in (|2.1(l yields (g, q, q, dq\j2^ M ^) as a 
local coordinate system on J 2 (M). We set q := dq\j2^ M y Thus, a local coor- 
dinate system q on M induces the natural local coordinate system (q, q, q) on 
J 2 (M). For further details on this affine bundle see ^Q, Q, [FT] . 

A Birkhoffian corresponding to the configuration manifold M is a smooth 
1-form u> on J 2 (M) such that, for any x G M, we have 

i> = (2.2) 
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where i x : f3^ 1 (x) — > J 2 (M) is the embedding of the submanifold /3 _1 (a;) into 
J 2 (M), (3 = ttm o ttj. From this definition it follows that, in the natural local 
coordinate system (q, q, q) of J 2 (M), a Birkhoffian u is given by 

m 

w = £0i(g,?,9)d^ (2.3) 
J=l 

with certain functions : J 2 (M) — > R. 

The pair (M, u>) is said to be a Birkhoff system (see 0). 

The differential system associated to a Birkhoffian to (see [7] ) is the set 

(maybe empty) D(u>), given by 

D(lu) := {z e J 2 (M) | oj(z) = 0} (2.4) 

The manifold M is the space of configurations of D(u), and D(u>) is said to have 
m 'degrees of freedom'. The Qi are the 'generalized external forces' associated 
to the local coordinate system (q). In the natural local coordinate system, D(w) 
is characterized by the following implicit system of second order ODE's 

Qj(q, q, q) = for all j = l,m (2.5) 

We conclude that the Birkhoffian formalism is a global formalism for the dy- 
namics of implicit systems of second order differential equations on a manifold. 
Let us now associate a vector field to a Birkhoffian u>. 

A vector field Y on the manifold TM is a smooth function Y : TM -> TTM 
such that Utm ° 5^=id. Any vector field Y on TM is called second order vector 
field on TM if and only if Tn M {Y v ) = v for all v € TM. 

A cross section X of the affine bundle (J 2 (M), ttj, TM), that is, a smooth 
function X : TM — » J 2 (M) such that 7rjoX=id, can be identified with a special 
vector field on TM, namely, the second order vector field on TM associated to 
X. Indeed, because (J 2 (M), ttj, TM) is a sub-bundle of (TTM, tttm, TM) as 
well as of (TTM,Tttm,TM), its sections can be regarded as sections of these 
two tangent bundles. Thus, using the canonical embedding i : J 2 (M) -> TTM, 
X can be identified with Y, that is, Y — i o X. 

In natural local coordinates, a second order vector field can be represented as 

m 

y = E 

A Birkhoffian vector field associated to a Birkhoffian u> of M (see 0) is 
a smooth second order vector field on TM, Y = ioX, with X : TAf J 2 (M), 
such that M C D(lu), that is 

X*w = (2.7) 

In the natural local coordinate system, a Birkhoffian vector field is given by the 
expression Ij2.6|l . such that Qj(q, q, q(q, q)) = 0. 



1 +1 W) ?)^t- 



(2.6) 
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A Birkhoffian u> is regular if and only if 

^0 (2.8) 

i,j=\,...,m 

for all (q, q, q), and for each (q, q), there exists q such that Qj(q, q, q) = 0, j = 
1, ...,m. 

If a Birkhoffian w of M is regular, then it satisfies i/ie principle of determinism, 
that is, there exists an unique Birkhoffian vector field Y = i o X associated to 
lu such that ImX = D(cu) (see [7]). 

A Birkhoffian to of the configuration space M is called conservative if and 
only if there exists a smooth function E u : TM — > R such that 

(X*w)Y = dE u (Y) (2.9) 
for all second order vector fields Y — i o X (see [J]). 

Equation l|2.9|l is equivalent, in the natural local coordinate system, to the 
identity (see 0, p. 16, eq. 4) 

m m 

J^Qjiq q, q)q ] = ^ 
i=i j=l 

is constant on TM if and only if dE UJ (Y) = for all second order vector 
fields Y on TM (see 0). 

If a; is conservative and Y is a Birkhoffian vector field, then 1)2. 9fl becomes 

dE w (y) = (2.11) 

This means that is constant along the trajectories of Y. 

We now introduce the concept of a dissipative Birkhoffian. 
A vertical 1-form on TM (see for example ^31) is a 1-form ^ on TM such that 
ty(V v ) = 0, for all V vector field on M, where V v is the vertival lift of the vector 
field V to TM. The local expression of a vertical 1-form is 

m 

* = $>iM)<V (2-12) 

3=1 

A 1-form D on TM is called dissipative if and only if D is vertical and D(Y) > 
0, for all Y second order vector field on TM. Allowing for (|2.12|) . the local 
expression of D is 

m 

D = Y / D J (q,q)dq 1 (2.13) 

3=1 

and from (|2.13|) . I|2.6|) . the inequality D(Y) > becomes 



det 



dqi 



(q, q, q) 



dEL 
dqi 



dE„ 



(2.10) 



m 



(2.14) 



A Birkhoffian lu of the configuration space M is called dissipative if and only 
if there exists a smooth function Eq^ : TM — » R such that 



(x*lo)y = «i£; 0u) (y) + D(y) 



(2.15) 



for all second order vector fields Y = i oX on TM, D being a dissipative 1-form 
on TM. 

Equation 12.15fl is equivalent, in a local coordinate system, to the identity 



3=1 3=1 



In view of (|2.14(l . we obtain from (|2.15|l . 

{X*lj)Y > dE 0ul (Y) 



(2.16) 



(2.17) 



for all second order vector field Y = i o X . That is equivalent, in local coordi- 
nates, to the dissipation inequality 



3=1 3 = 1 



1 H — ^-—9 



9gJ 



9g3 



(2.18) 



If a; is a dissipative Birkhoffian and Y is the Birkhoffian vector field, then l|2.17|l 
becomes 

dE „(Y)<0. (2.19) 

This means that Eq^ is nonincreasing along the trajectories of Y. 
If the dissipative 1-form on TM has the particular expression 



i,3=l 



(2.20) 



when we calculate the function D(Y) on TM, we obtain the so called Rayleigh 
dissipation function 1Z : TM — > R 



(2.21) 



3 RLC circuit dynamics 

A simple electrical circuit provides us with an oriented connected graph, that is, 
a collection of points, called nodes, and a set of connecting lines or arcs, called 
branches, such that in each branch is given a direction and there is at least 
one path between any two nodes. A path is a sequence of branches such that 
the origin of the next branch coincides with the end of the previous one. The 
graph will be assumed to be planar, that is, it can be drawn in a plane without 
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branches crossing. For the graph theoretic terminology, see, for example 
Let b be the total number of branches in the graph, n be one less than the 
number of nodes and m be the cardinality of a selection of loops that cover the 
whole graph. Here, a loop is a path such that the first and last node coincide 
and that does not use the same branch more than once. By Euler's polyhedron 
formula, b — m + n. 

We choose a reference node and a current direction in each /-branch of the 
graph, I — 1, b. We also consider a covering of the graph with m loops, and 
a current direction in each j-loop, j — 1, m. We assume that the associated 
graph has at least one loop, meaning that m > 0. 

A graph can be described by matrices: a (6n)-matrix B € 3Jt(,„(R), rank(B) = 
n, called incidence matrix and a (6m)-matrix A G 9Ttf, TO (R), rank(yl) = m, 
called loop matrix. These matrices contain only 0,1, -1. An element of the 
matrix B is if a branch b is not incident with a node n, 1 if branch b enters 
node n and -1 if branch b leaves node n, respectively. An element of the matrix 
A is if a branch b does not belong to a loop m, 1 if branch b belongs to 
loop m and their directions agree and -1 if branch b belongs to loop m and their 
directions oppose, respectively. For the fundamentals of electrical circuit theory, 
see, for example [5]. 

The states of the circuit have two components, the currents through the branches, 
denoted by I e R b , and the voltages across the branches, denoted by v € R h . 
Using the matrices A and B, KirchhofFs current law and Kirchhoff's voltage 
law can be expressed by the equations 

B T \ = (KCL) (3.1) 

A T v = {KVL) (3.2) 

Tellegen's theorem establishes a relation between the matrices A T and B T : the 
kernel of the matrix B T is orthogonal to the kernel of the matrix A T (see e.g., 
4, page 5). 

The next step is to introduce the branch elements in a simple electrical 
circuit. The branches of the graph associated to a RLC electrical circuit, can 
be classified into three categories: resistive branches, inductor branches, and 
capacitor branches. Let r denote the number of resistive branches, k the number 
of inductor branches and p the number of capacitor branches, respectively. We 
assume that just one electrical device is associated to each branch, then, we have 
b = r + k + p. Thus, we can write I = (ipj, i( a ), i Q ) 6 R r x R fe x R p ~ R b , where 
Irpi, i( a ), i a are the currents through the resistors, the inductors, the capacitors, 
respectively, and v = (up] , Vt a ), v a ) G R r x R fe x W ~ R b , where vm, V( a ), 
v a describe the voltage drops across the resistors, the inductors, the capacitors, 
respectively. 

Each capacitor is supposed to be charge-controlled. For the nonlinear capacitors 
we assume 

t'o = Ca(Q a ), a=l,...,p (3.3) 

where the functions C a : R — ► R\{0} are smooth and invertible, and the Q a 's 
denote the charges of the capacitors. The current through a capacitor is given 
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by the time-derivative of the corresponding charge 



la = ~rjf> a = l,...,p (3.4) 

t being the time variable. 

Each inductor is supposed to be current-controlled. For the nonlinear inductors 
we assume 

v a = L a (i Q )-A a = l,...,k (3.5) 
at 

where L a : R — ► R\{0} are smooth invertible functions. 

There are several types of nonlinear resistors, among them current controlled 
resistors and voltage controlled resistors. Generally, their constitutive relations 
are defined by some continuous functions of Ir and vr , that is, 
/r(i r ,fr) = 0,r = l,...,r. 



3.1 Current controlled resistors 

We first consider the case that the nonlinear resistors are current controlled, 
that is, the constitutive relations are given by 



vr=Rr{lr), T = l,. 



(3.6) 



where Rr : R — ► R are smooth functions. In order to obtain a dissipative 
Birkhoffian, we assume that, for all i^O, 



Rr{x)x>0, Vr = l,...,r 



(3.7) 



that is, for each nonlinear resistor, the graph of the function Rr lies in the union 
of the first and the third quadrant. 

For linear resistors the relations (|3.6|> can be written in the form 



vr = Rrir 



(3.8) 



where Rp > are real constants. 

Taking into account l|3.3ll . (|3.4I) . 13.5(1 . I|3.6[) . the equations (|3.1|l . (|3.2|l become 



B 1 




i?r(ir) 
C a {Q a ) 



(3.9) 



In the following we give a Birkhoffian formulation for the network described by 
the system of equations ]3.9\) . using the same procedure as in '81. That is, using 
the first set of equations IH.9\) . we are going to define a family of m- dimensional 
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affine-linear configuration spaces M c C R & parameterized by a constant vector 

c in R". This vector is related to the initial values of the Q-variables at some 

instant of time. A Birkhoffian lo c on the configuration space M c arises from a 

linear combination of the second set of equations y8.y\) . Thus, (M c ,u c ) will be 

a family of Birkhoff systems that describe the considered RLC network. 

We notice that the first set of equations (|3.9|) remains exactly the same for linear 

and nonlinear electrical devices. Thus, for obtaining the configuration space, it 

is not important whether the devices are linear or nonlinear. 

Let H : R b — > R" be the linear map that, with respect to a coordinate system 

(x 1 , x b ) on R b , is given by 



H{x l ,...,x h ) = B 



T 



( * 



\ 



(3.10) 



Then, 7J -1 (c), with c a constant vector in R™, is an affine-linear subspace in 
R b . Its dimension is m = b — n, because rank(B) = n. 
We define M c as 

M C :=H- X (c) (3.11) 

We denote local coordinates on M c by q = (q 1 , .., q m ). Then, the natural coor- 
dinate system on the 2-jets bundle J 2 (M C ) is given by (q,q,ij). 
We will now represent the Birkhoffian in a specific coordinate system on M c : 
In the vector spaces R r , R fc , we identify points and vectors 

dQ[r] dQ {a) 

Ir -=^T' Ia -=^r (3 - 12) 

with (Q[r])r=i,...,rj (Q( a ))o=i,...,fc coordinate systems on R r , respectively, on R fc . 
Taking into account (|3.12|l and the fact that the matrix B T is a constant matrix, 
we integrate the first set of equations (|3.9I) to arrive at 

, Q [r] 

B T Q (a) I =c (3.13) 



with c a constant vector in R". 

Likewise consider coordinates in R b ~ R r ^ x R fc x R p , 

x :— Q[ X ], x r := Q[ r ], x r+ := Q (1) , x r+ := Q (fe) , 

:= Qii ^ x b ._ Qp ( 314 ) 

From H3.10(l . H3.11|l . we see that we can define coordinates on M c by solving 
the equations 13.13|l in terms of an appropriate set of m of the Q-variables, say 
q = (q 1 , q m ). In other words, we express any of the x-variables as a function 
of q = (q 1 , ...,<7 m ), namely, m 

x r = y^^J\ffqf + const, T = 1, r 

3=1 
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x a = y^A/"jV + const, a = r+l, ...,r + k 

m 

a; Q =^JV^ g J ■ + const, a = r + k + l,...,b (3.15) 

i=i 

with certain constants AIJ, J\ff, Afj*. Here we can think of the constants const 
as being initial values of the ^-variables at some instant of time. 
From (EH, (|3~TH . lj3~Hjl and differentiating (|3~15|) we get 

l = Mq (3.16) 

with the matrix of constants N G 9Jtb m (R), for some q G R m . 

Using Tellegen's theorem and a fundamental theorem of linear algebra, we now 

find a relation between the matrices N and A. 

By a fundamental theorem of linear algebra we have 

[Ker{A T ))- L = Im{A) (3.17) 

where A G 2tt 6m (R), A'er(A T ) := {x G R b | A T x = 0} is the kernel of A T , 
/m(yl) := {x G R b | Ay — x, for somey G R m } is the image of A and - 1 denotes 
the orthogonal complement in R b of the respective vector subspace. 
For the incidence matrix B G £E7t&„(R) and the loop matrix A G 9Jtfc m (R), which 
satisfy Kirchoff's law l|3.1ll . (|3.2(l . Tellegen's theorem writes as 

Ker(B T ) = {Ker{A T )) ± - (3.18) 

From the first set of equations in H3.9|l . and by constraction of the matrix J\f in 
H3.16|l . we have 

Ker(B T ) = Im(N) (3.19) 

Therefore, using (|3~T7|l . l|3~TH|l . lj3~TU|l . we obtain 7m(A) = Im(Af). Then, 
another application of l|3.17|l yields 

Ker(A T ) = KeriAf 7 ) (3.20) 

Taking into account (|3.20|l . we see that there exists a nonsingular matrix £ G 
Wt mm (R) satisfing 

€A T — J\f T (3.21) 

The matrix £ provides a relation between the vector of the m independent loop 
currents and the coordinate vector q introduced on M c . 

Taking into account \3.2U\) . we define the Birkhoffian uo c of M c such that the 
differential system y2.5}) is the linear combination of the second set of equations 
in l^y.yp obtained by replacing A T with the matrix M T . Thus, in terms of q- 
coordinates as chosen before, the expressions of the components Qj(q, q, q) are 

Q j (q,q,q) = F j (q)q + H j (q) + G j (q), j = l,...,m (3.22) 



10 



wher 



Fjm-= E ^.-r ^ =E E ^Vr(?)? ( 3 - 23 ) 

a = r + l \ Z= 1 / \i = l / £=1 a = r + l 

r / m \ ?' 

^ (,) := e <«r EX*' =E w) ( 3 - 24 ) 



E N?C a _ r ^ k (q) (3.25) 



G 3 -( g ):= E *?C a - r - k E^V 



const 



We note that the Birkhoffian 1)3.22(1 is not conservative. We easily see that 
there does not exist a function E u such that 1)2.10)1 is fulfilled for the Birkhoffian 

E23> since w ^ w- 

For an i?LC electrical network with nonlinear resistors, described by y3.b}) . 
\3. 7| ), nonlinear inductors and capacitors described by $3.5\) . respectively \3. 3\) . 
we claim that the Birkhoffian 13. 221) is a dissipative Birkhoffian. 

Indeed, in the view of the assumption 1)3. 7J) . the vertical 1-form Z? on TM 
given by 



3=1 

with Hj(q) in 13.24|l . is dissipative, that is, 



3=1 



r=i 



(3.26) 



(3.27) 



We showed in jS] that the following smooth function Eq w on TM 



r + k m m 

E E E 

a— r + 1 l — l ii<...<ii—l 



{l-l)L(t 2) {q)]N? 1 ...N? l dq i K..d4 i < + 



E E E / c^ZrUiiW^-K^ 11 -^ 11 ( 3 - 28 ) 

a=r+fc+l !=1 ii<...<ii = l w_/ 



satisfies the identity 



^—9 + -^rri 



(3.29) 



According to (|3.29|l . the Birkhoffian l|3.22|l satisfies 1)2.16)1 with the function 
E 0ui (q, q) given by H3.28fl and the dissipative 1-form D given by (|3.2(i|) . ■ 

Let us now discuss the regularity of the Birkhoffian given by 13.2211 . 

If there exists in the network at least one loop that contains only capacitors, 
or only resistors, or only resistors and capacitors, then the Birkhoffian 
associated to the network is never regular. 
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In [5] we have shown that if there exists at least one loop in an LC network 
that contains only capacitors, then the Birkhofnan associated to the network is 
never regular. The Birkhofnan associated to an RLC network which contains 
at least one loop formed only by resistors or only by resistors and capacitors is 
never regular as well. The proof is based on the fact that for the Woop which 
does not contain any inductor branches, on the column I of the matrix A we 
have 

Af = 0, a = r + l,...,r + k (3.30) 
For the Birkhofnan i|3.22[l . the determinant in (|2.8|l becomes 



det 



dQj 

dif 



det 



i.j — l,...,m 



r+k 



a— r+1 



(3.31) 



- i,j — l.....m 



From (|3~2T|) . we get Af? = YX=i A h > for any a = r + 1, r + k and taking 
into account (|3.30() . we have 



r+k 



i 1= i 



r+k 



]T (Al) 2 L a _ r (q) 



_a=r-\-l 



£ (q ef + et 4 



i l<31 
<1 ,317*1 



r+k 



J2 AlAlL a . r (q) 

(3.32) 



_a— r+1 



Using basic calculus, the determinant of the matrix with the elements (|3.32f) . 
is a linear combination of determinants having at least two linearly dependent 
columns. This shows that the determinant in the right hand side of (|3.31(l is 
equal to zero. Thus, the Birkhofnan H3.22(l is not regular. ■ 

We now discuss the question, how to proceed in the case than the Birkhofnan 
given by l|3.22|l is not regular in the sense of definition (|2.8|) . 

// there exists in the network m\ < m loops which contain only capaci- 
tors, all the other loops containing at least an inductor, we can regularize the 
Birkhoffian yj.2fy) via reduction of the configuration space. The reduced 
configuration space M c of dimension m — m,\ 7 is a linear or a nonlinear subspace 
of M c , depending on whether the capacitors are linear or nonlinear. We claim 
that the Birkhoffian uj c on the reduced configuration space M c is still a dissipa- 
tive Birkhoffian. Under certain conditions on the functions L a , a = 1, ...,k, 
which characterize the inductors, the reduced Birkhoffian Co c will be a regular 
Birkhoffian. 

Without loss of generality, we can assume that there is one loop in the network 
that contains only capacitors and in the coordinate system we have chosen 



Af[ = o, r = i, . 



TV? = 0, a = r + l, 



(3.33) 



Thus, the Birkhoffian components il3~25|> . with (|3~21) . (|3~25|) . are given 

by, j = 2,...,m, 
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b 

Ql(<Z, <?,<?) = E X?C a _ r _ k {q) 

a=r-\-k-\-l 

in r+fc r b 

Q 3 (q,q,q) = E E K^ 1 — (?) ? + E (4) + £ ■A/J , 6 01 _ r _ fc (g) 

i=2o=r+l T=l a=r+fe + l 

(3.34) 

We note that, according to l|3.33[l . q 1 does not appear in any function Rr(q), 
L a -r(q) and the terms L a _ r (g)g 1 do not appear in any of the BirkhofRan com- 
ponents Q 2 (q, q,q),--, Q m (q, q,q). 

We define the (m — l)-dimensional nonlinear space M c C M c by 

b 

M c ={qeM c \ J2 ^iC a - r -k(q) = 0} (3.35) 

a=r+fc+l 

By the implicit function theorem, we obtain a local coordinate system on the re- 
duced configuration space M c . Taking q 1 :— q 2 ,..., q" 1 ^ 1 :— q m , the Birkhoffian 
has the form uj c = X^Jli 1 Qjdq^ , 

Qj(q, ?. £~) = Fj(q)q + Hj(q) + Gj(q), where (3.36) 

m — 1 r+fc /m—1 \ 



urn ■■= E E ^s+D^+D^-r E ^?+Dff ( 3 - 37 ) 

i=l a =r+l \ i=l / 

r /m—1 \ 

:= E^O+D^ E ^+1)9 ( 3 - 38 ) 
r=i V (=i / 

b I m-1 \ 

c,(g) := E Mifti 1 ' •••> 9 m_1 ) + E ^P+i)^ + 

a=r+fc+l \ 1=1 J 

(3.39) 

/ : U C R" 1 ^ 1 — > R being the unique function such that /((fa) = <7ch <7o ^ R, 
and 

6 / m-l \ 

E M?C a -r-k Kf(q\ ■■; 3™^) + E "^(W + co ^ = ( 3 - 40 ) 
u=r+k+l V i=l / 

for all q = (q 1 , q" l ~ x ) S C/, with IT a neighborhood of go = {qh, So* - )• 
We will now prove that the Birkhoffian (|3.36(l is dissipative. In order to do so, 
we will show that there exists a function Eq^ (g, q) satisfying 



E 9, = E 

3=1 3=1 



9Eq „ j = ,_ „ 



(3.41) 



where I? = X^i 1 ^ifo ( ?)^'?' ; i s a dissipative l-form on TM C 
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We consider the following 1-form on TM C 

rn — 1 r 



j=i r=i 



(3.42) 



In the view of the assumption l|3.7|l , the vertical 1-form 13.42|l on TM C is dissi- 
pative, that is, 



£ 

r=i 



m — 1 



'm-l 

E^W 



> 



(3.43) 



Therefore, (|3.41|) is fulfilled if Eq^ (q, q) can be chosen in such a way that 

m — 1 m — 1 



d£ , ._, | ago, 

9gJ 9 dq J 



(3.44) 



is satisfied. Because of the special form of the terms on the left side of 
we may assume that Eq^ (q, q) is a sum of a function depending only on q, and a 
function depending only on q. From the theory of total differentials, a necessary 
condition for the existence of such functions is the fulfillment of the following 
relations 

df, (g) _ dfi(q) _ q 



dq aq 3 
dGjjg) _ dOgg) = 

dq 1 dq J 



(3.45) 



for any j, I = 1, .. . , m — 1 , where 

m— 1 r+A; 



f m — 1 



^■(9) : = E E ^ +1) ^ +1) L a _ r £^? +1) ? 



(3.46) 



i— 1 a— r+1 

From H3.46JI . we get: 

r+k 



1 = 1 



dm) 



rn—l r-\-k 



#9 a=r+l i=l a=r+l 

(3.47) 

where L' a _ r := dLa ^Z^ • Then, the left side of the first relation in H3.45J) becomes 

r+fe 

E [K+d^+d -^3+1)^5+1)) 

a— r+1 

- -^+1)^+1)) ( E ^Ti+l)^) Z a-r(€ (3-4 
i=l 
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We easily see that the expression in (|3.48[) is zero, thus the first relation in l|3.45[l 
is fulfilled. 

From i|3.39[l . the second relation in (|3.45|1 reads as 



E \^A-r-M 
a=r+k+l 

AA(7 +1 )C„_ r _ fe (g) 



= 



(3.49) 



where C' 



r _ k : = dc a -r- k (v) _ xhe re i a tion ijS^j) reduces to 



a=r+/c+l 



(3.50) 



Taking into account (|3.4U|) . the above relation is fulfilled, for any j, I — 1, to— 
1. Indeed, taking the derivatives with respect to cp and also to q l , in the equation 
H3.40f> . we obtain, respectively, 



ot=r-\-k-\-l 
b 

E ^?^-r-fc(«) 



7Vl +yv 0'+i) 



a+i) 



= o 



= o 



(3.51) 



Now we multiply in l|3.51|l the first equation with 9 g-f , the second equation 
with — d gqP and we add the resulting equations to obtain the equation (|3.5Q(I . 

Thus, we proved the existence of a function E 0aj (q,q) such that (|3.41l) is 
fulfilled, with the dissipative 1-form given by (|3.42(l . Therefore, the Birkhoffian 
H3.36JI is dissipative. 

For the Birkhoffian Ij3.36|l . the determinant in (|2.8|l becomes 



det 



(9, k, q) 



det 



— l,...,m — 1 



,a=r-\-l 



(3.52) 



is 



If the determinant in (|3.52|) is different from zero, then the Birkhoffian 
regular. ■ 

// there exists in the network 7712 < m loops which contain only resistors, all 
the other loops containing at least an inductor, we can regularize the Birkhof- 
fian IS. 22\) via reduction of the configuration space. The reduced config- 
uration space M c of dimension m — rri2, is a linear or a nonlinear subspace 
of M c , depending on whether the resistors are linear or nonlinear. We claim 
that the Birkhoffian uj c on the reduced configuration space M c is still a dissipa- 
tive Birkhoffian. Under certain conditions on the functions L a , a = l,...,fc, 
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which characterize the inductors, the reduced Birkhoffian uj c will be a regular 
Birkhoffian. 

Without loss of generality, we may assume that we have one loop in the 
network that contains only resistors and in the coordinate system that we have 
chosen, the constants read as 

A/"" = 0, a = r + 1, r + k, = 0, a = r + k + 1, ...,b (3.53) 

I) Let us first consider the case that the resistors in this loop are linear 
resistors, that is, described by (|3.8|) . all the other electrical devices in the network 
being nonlinear. This means that we have 

A/fVO, T = l,...,r im , Af[ = 0, T = run + 1, r (3.54) 

where ru n is the number of linear resistors in the network. 

In this case, the expressions of the Birkhoffian components i|3.22[) . with (|3.23|) . 

(|3.24[) . (|3.25|) . are given by, j = 2, ...,m, 



r=i 1=1 

rn r-\-k T'Un m 

i=2 a=r+l r=l 1=1 

r b 

£ ^Rr{q)+ £ NfC a - r - k {q) (3.55) 

T=rn n +1 a =r+k+l 

We note that according to l|3.53|l . (|3.54|) . q 1 does not appear in any function 
Rriq), L a - r {q), the terms L a _ r (g)g 1 do not appear in any of the Birkhoffian 
components Q2(q, q, q), Qm(q, q, q) as well q 1 does not appear in any function 

C a - r -k{q)- 

We define the (m — l)-dimensional linear space M c C M c by 

run m 

M c = {q e M c I £ J2^i R ^1 l + c i = °} ( 3 - 56 ) 

r=i i=i 

with ci a real constant. 

We take q 1 := q 2 ,---, q" 1 ^ 1 := q m , as local coordinates on the reduced configu- 
ration space M c . Then, making use of [|3.56JI and of the fact that M[ ^ and 
Rr > 0, for any T — l,...,ru n , we can express q 1 as a linear combination of 
q , ...,q , denoted g(q), such that 



r=i 



m— 1 

Mi 9(4) + E^G+i)? 



i=i 



= (3.57) 



Thus, the reduced Birkhoffian has the form ui c = X)j=i Qjdq 3 , 

Qj{q, q, 'q) = Fj(4)'q' + Hj(q) + Gj(q), where (3.58) 



1G 



tn—1 r-\-k 



i— 1 a— r+1 



(3.59) 



r=i 



JVf <?(<?) + E ^+1)9 



z=i 



f m — l 



E ^ +1) i? r (E^ ( W 

r=n,„ + l 



. 1=1 



Gj(q) 



b (m-1 \ 

£ 7V5 +1) c a _ r _ fc E ^a+D^ + const 

a=r+k+l \ 1 = 1 I 



(3.60) 



(3.61) 



The BirkhofRan given by H3.58|l is still dissipative. We will see that there exists 
a function Eq u (q, q) such that 



m — l 



E & <?)^ = E 

3=1 3=1 L 



9 H — <7 +Dj(q,q)q 



dq3 



dq 



where £> = X^i 1 Dj{q, q)dq^ is a dissipative 1-form on TM C . 
We consider the following 1-form on TM C 



(3.62) 



D=Y,H j (q)dq i 

3 = 1 

Let us check that the vertical 1-form l|3.63|l is dissipative, that is, 

m — 1 

Y^Hjm 3 >o 

3=1 

From H3.60fl , the left side of l|3.64|l writes as the sum S\ + S2 , where 



m—l run 

3=1 r=i 



jVlg(q) + £ ^(i+i) 



S 2 = e 

r=r H „+i 



( E J^+i) 

3 = 1 



3=1 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



We now multiply the equation l|3.57|l by the function g(q). Using the resulting 
equation we can write the sum in (|3.65|) in the form 



Si = Y R r 



r=i 



n 2 



Ml 9(4) + E M( l+ i } q 



1=1 



(3.67) 
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Since Rr > 0, T = 1, ru n , the sum Si is strictly positive. 

Because all nonlinear resistors considered satisfy the condition l|3.7[l . the sum S2 

in (|3.66|) is strictly positive as well. Therefore, the inequality (|3.64|) is fulfilled 

and the vertical 1-form in i|3.63[) is dissipative. 

We now look for a function Eq^ (q, q) such that 



m — 1 



X7 

q 



3=1 



E 

3=1 



9E -j 



(3.68) 



Because of the special form of the terms on the left side of (13.6811 . we can look 
for the required function Eq^ (g, q) as a sum of a function only depending on <J, 
and a function only depending on q. From the theory of total differentials, a 
necessary condition for the existence of such functions is the fulfilment of the 
following relations 



(4) „ dFijq) = q 
dq d'q 1 



(3.69) 



dG^q) _ dG,(q) _ 



dq 1 



dqi 



for any j, I = 1, m — 1, where 



tn—X r-\-k /m—1 \ 

m) = EE ^ +1) ^ +1) L Q _ r "fad i 

i=l a=r+l \ 1=1 J 



(3.70) 



From (f^TTOfl . (|3~HT|) we get 



dq 



r+A: 



a— r+1 



m — 1 r+A; 

E E 

2=1 a— r+1 



E ^( j+ l)^(l+l)La-r (€ + E E ^5+1)^+1)^+1}^^)^ 



ggj(g) 
dq 1 



— E ^§+l)^?+l)^a-r-fe(#) 



(3.71) 
(3.72) 



a— r+A;+l 



where L' 



dL a - r (ri) ^1 

dr\ ' a— r— ft 

are fulfilled. 



d?7 



We can easily check that the 



a—r 

equations 

Thus, we proved the existence of a function Eq^ (q, q) such that l|3.62[l is fulfilled, 
with the dissipative 1-form given by (|3.63(l . that is, the Birkhofnan 1)3.58(1 is 
dissipative. 

For the Birkhofnan (|3.58|l . the determinant in 1)2. 8|l becomes 



det 



dQj 



(<7, q, q) 



det 



— l,...,m — 1 



r+fc 



a— r+1 



(3.73) 



If the determinant in (|3.73() is different from zero, then the Birkhoffian (|3.58() is 
regular. ■ 
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II) Let us now consider the case that the resistors in the loop formed only 
by resistors are nonlinear devices too, that is, they are described by l|3.6[) . with 
the assumption i|3.7fl . Now, instead of i|3.54|) we have Af[ ^ 0, T = 1, r. 
The component Qi(q, q, q) of the BirkhofHan takes the form 

r 

Ql(q,q,q) = J2^lMq) (3.74) 
r=i 

and the other components Q2(q,q,q),---, Qm(q,q,q) are the same as in (|3.55H 
with the terms following ru n = absent. According to (|3.53|) . q 1 does not ap- 
pear in any function L a - r (q), the terms La-^qjq 1 do not appear in any of the 
BirkhofHan components Q2(q, q,q), Qm(q, q, q) as well q 1 does not appear in 
any function C a - r -k{q)- 

Using (|3.74|) . we intend to define the (m — l)-dimensional configuration space 
M c C M c . The relation l|3.74[l is a nonlinear velocity constraint, which in general 
is a nonholonomic constraint. Nevertheless, because of this constraint imposed 
on the system, the equations which describe the dynamics are 
Q2(q, q, q) = 0,..., Q m (q, q, q) = 0. 

Taking q 1 :— q 2 ,..., q" 1 ^ 1 := q m , a coordinate system on the reduced configura- 
tion space M c , the Birkhoffian has in this case the form Cj c = X}j=i Qj^Fj 

m — 1 r+k /m—1 \ 

QMll) = E E ^U + l)^ + l)La-r E ^8+1)9 K + 
i=l a =r+l \ 1=1 / 

r / m — 1 \ 

E^ + D^r U?h(i\..J m )+ E^ + x)9 
r=i V i=i I 

b fm-1 \ 

+ E Mu+i)C a -r- k E^+i)«M ( 3 - 75 ) 

a=r+k+l \ 1=1 / 

where h : U C R" 1 ^ 1 — ► R is the unique function such that h(q Q ) — <Jq, q\ e R, 
and 

r / m—1 \ 

E<i? r k^^rVEW =° (3J6) 

r=i \ i=i ) 

• • 1 • Tfi — 1 • * 1 • Tfl — 1 

for all q = (q , q ) £ [/, with U a neighborhood of q = (q , % ). 
One can prove using the same ideas as in the previous case, that the Birkhoffian 
given by l|3.75[) is still dissipative. 

If the determinant det Z)oir+i ■ /v (j+i)- ,v (i+i) J 
the Birkhoffian (|3.75|) is regular. 

// there exists in the network < m loops which contain only resistors and 
capacitors, all the other loops containing at least an inductor, we can regularize 
the Birkhoffian yd. 2ty) by introducing into each of these 7713 loops, an inductor in 
series, with the inductance functions £ a > ■ R — > R\{0}, a' = 1, ...,7713, having 



E::; +1 ^. +1) 4 +1) L _ r M + 0, then, 

\ / J i,j=X,,,,,rn— 1 
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very small values. The configuration space remains M c of dimension m. under 
certain conditions on the functions C a > and on the functions L a , a = 1, ...,k, 
which characterize the others inductors, the Birkhoffian io ext on M c will be a 
dissipative regular Birkhoffian. 

Without loss of generality, we may assume that we have one loop in the 
network containing only resistors and capacitors and in the chosen coordinate 
system 

A/*" = 0, a = r+l,...,r + k (3.77) 
The expressions (|3.22|) of the Birkhoffian components become, j = 2, ...,m., 
Qi(q,q,q) = j2M?Rr(q)+ ]T JV/?C a _ r _ fc (g) 

r=l a=r+fc+l 
m r-\-k r b 

Qj(q, q, q) = £ £ ^iK^a-r (?) ? + J2 M F^ (?) + E 

j=2 a=r+l T=l a=r+k+l 

(3.78) 

We introduce into this loop an inductor in series, described by the following 
relation between the current and the voltage 

v=d(i)j t (3.79) 

L\ : R — > R\{0} being smooth invertible function. After introducing the 
inductor, the number of branches of the graph associated to the circuit increases 
by one, that is, there will be b + 1 branches, and the number of nodes increases 
by one as well, that is, n becomes n + 1. Still the cardinality of a selection of 
loops which cover the whole graph remains m. The configuration space is the 
same M c , with dimension m. The corresponding Birkhoffian, denoted by oj^ xt , 
has the component Qi(q, q, q) given by, 

r b 

Q 1 (q 7 q,q)^C(q 1 )q 1 + J2^Mq)+ ]T JV?C„_ P _ fc (g) (3.80) 

T=l a=r+k+l 

the others Q2(q,q,q),..-, Q m (q,q,q) have the same form as in (|3.78|) . This 
Birkhoffian is dissipative, the expression of the function E 0ui on TM C is (|3.28|l 
plus the term J C(q 1 )q 1 dq 1 . The dissipative 1-form has the form (|3.26l) . 
If Ciq^det [E^+i La-r (q)} . . , , ¥= 0, this Birkhoffian is 

regular. ■ 



i 7 j—l,...,m—l 



3.2 Voltage controlled resistors 

Let us now consider the nonlinear resistors for which the constitutive relations 
are given by 

ir = £Kr(wr)» r = l,...,r (3.81) 
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where 9tr : R — > R are smooth functions. In order to obtain a dissipative 
Birkhoffian, we also assume that, for all x ^ 0, 

m r (x)x>0, Vr = l,...,r (3.82) 

that is, for each nonlinear resistor, the graph of the function 9V lies in the union 
of the first and the third quadrant. 

Taking into account (|3.3|l . (|3.4|l . I|3.5|l . I|3.81|l . the equations (|3.1|l . (|3.2[) govern- 
ing the circuit have the form 




< f / vr \ (3.83) 

A T L a {l a )ia -0 
< V C a (Q a ) J 

[ ir^^rM 

Using the first set of equations \S.H!^) . we define a family of m- dimensional 
affine-linear configuration spaces M c C R b parameterized by a constant vector 
c in R™. This vector is related to the initial values of the Q-variables at some 
instant of time. A Birkhoffian to c on the configuration space M c arises from 
the second set of equations \3.83\) . Thus, (M c ,lu c ) will be a family of Birkhoff 
systems that describe the considered RLC network. 

The first set of equations l|3.83|) is the same as the first set of equations H3.9(l . 
Thus, using them we can define the family M c of m-dimensional affine-linear 
configuration spaces (|3.11|) . For a coordinate system q = (q 1 , q m ) on M c , 
the relations between the x-coordinates (|3.14|) and the g-coordinates are given 
by H3.15fl . The matrix of constants TV satisfies (I3.16|) . Taking into account 
H3.20[l , we define the Birkhoffian lo c of M c such that the differential system (|2.5Jl 
is the linear combination of the second set of equations in (|3.83(l obtained by 
replacing A T with the matrix Af T . In terms of the g-coordinates chosen before, 
the components Qj(q,q,q) of the Birkhoffian have the implicit form 

!Qj(q, Q, <?) = Fj^q + Gjiq) +Er=i- /V 7 w r, j = l,...,m 
(3.84) 
Er=i^ i= ^r(«r) 

where the functions Fj(q)q, Gj(q) are given by H3.23|l . (|3.25|l . 
The Birkhoffian \3.8$ is dissipative . 
Indeed, with the function Eo uj (q,q) given by (|3.28l) . the identity l|2.16|) be- 
comes for the Birkhoffian (|3.84() . 

rn m 

j=i i=i 



^NjvTq* 



(3.85) 



r=i 
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It remains to show that the vertical 1-form defined implicitly by 



D j(<l> 4) = Er=i-^7 u r, j = l,...,m 
is dissipative. From the second set of relations in H3.86jl . we have 

m r r 

N J V ^ = 12 Wr(vr)v r 
j=i r=i r=i 



(3.86) 



(3.87) 



Therefore, the inequality l|2.14|l reads as 2r=i ^r(vr)vr > 0. The last inequal- 
ity is satisfied in view of the assumption l|3.82[l . ■ 

As in the case of current controlled sources, the Birkhoffian \3.84\) is never 
regular if the network contains closed loops formed only by capacitors, or re- 
sistors, or both of them. 



4 Examples 

Example 1): This example is based on the following oriented connected graph 




3 Figure 1 



We have r = 1, fe = 3, p = 3, n = 3, m = 4 , b = 7. We choose the reference 
node to be V4 and the current directions as indicated in Figure 1. We cover the 
associated graph with the loops h, h, h, h- Let V = (Vi, V2, V3) £ R 3 be the 
vector of node voltage values, 1 = (ijr^ifa)^) G R 1 x R 3 x R 3 be the vector 
of branch current values and v = {vrpuVra), v a ) e R 1 x R 3 x R 3 be the vector 
of branch voltage values. 

The branches in Figure 1 are labelled as follows: the first branch is the resis- 
tive branch Ri, the second, the third and the fourth branch are the inductive 
branches Li, L2, L3 and the last three branches are the capacitor branches Ci, 
C2, C3. The incidence and loop matrices, B G 9Jl73(R) and A e £Dl74(R), write 
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as 



B 



( 


-1 





1 


\ 




f 











-1 \ 










-1 














1 


-1 







-1 














1 


-1 










-1 


1 




, A = 










1 










1 











-1 


1 










1 














1 








-1 


\ 


-1 


1 





/ 




V 


1 








o / 



(4.1) 



One has rank(_B) = 3, rank(A) = 4. 

All the electrical devices are considered to be nonlinear and described by the 
relations, {£3, with Ci,C 2 ,C 3 : R -> R\{0}, Lx,L 2 ,L 3 : R -> R\{0} 

smooth invertible functions and by the relation (|3.fi|l , with R\ : R — > R smooth 
function, such that (j3.7(l is satisfied, that is, for any i / 



Ri(x)x > 

The equations l|3.9|l which govern the network have the form 



(4.2) 



-i[i] - 

-1(2) 

im - i 



Q 2 - 
-1(3) 



Q 3 = 

f Qi + Q 3 = 



[1] - 1(1) + 1(3) 



= 



-C 1 (q 1 ) + C 2 (q 2 ) + C 3 (q 3 ) = 

^2(l(2))i(2) +Ci(Q!) =0 

^l(l(l))i(l) - ^2(l(2))i(2) + ^3(l(3))i(3) 

-i?i(i[i]) - ii(i(i))i ( i) - C 2 (q 2 ) = 



The relations H3.12J1 . I|3.14(l read as follows for this example 



(4.3) 



Q[i]> 



i[i] := Q[i] 

3 



i( Q ) := 



Q(i 



Q(2): 



1(a), 



Q(3)) 



a = 1,2,3 



Qi, x 



Q 2 , x 



(4.4) 

Q 3 
(4.5) 

Using the first 3 equations of the system l|4.3l) we define the 4-dimensional 
affine- linear configuration space M c . In view of the notations l|4.4l) . (|4.5I) . we 
integrate these 3 equations and solving them in terms of 4 variables, we obtain, 
for example, x 2 — x 1 + x 4 + const , x 5 — x 3 +x 4 — x 7 + const, x 6 — x 1 +x 7 +const. 
Thus, a coordinate system on M c is given by 

„i 



7 2 

x',q 



x 4 , q 3 :— x , q 4 



(4.6) 



23 



The matrix of constants M = \ Nj 



from H3.15fl is given by 



r = l. 2,0 = 3,4,0 = 5,6,7 
3 = 1,2,3,4 



Therefore, in terms of the g-coordinates 1)4.6(1 . we 



/ 1 \ 
110 
1 

/V = oioo 

-110 1 
10 10 
\ 1 / 

may define the Birkhoffian oj c = Qi(g, q, q^dq 1 + Q 2 (q, g, q)dq 2 + Qs(g, q, q)dq 3 + 
Q 4 (g,9,g)dg 4 of M c as in with QOgl , EH5| l. that is, 

Qi(g, g, g) = -C^-g 1 + g 2 + g 4 + const) + C 2 (q 1 + g 3 + const) + C^q 1 ) 

Q%{q, q, q) = [Li(q 2 + f) + L 3 (q 2 )] g 2 + L : (g 2 + g 3 )g' 3 + C^-g 1 + g 2 + g 4 + const) 

Q 3 (q, g, g) = Li(g 2 + g 3 )g 2 + Li(g 2 + g 3 )g 3 + ^(g 3 ) + C 2 (q 1 + g 3 + const) 

Q 4 (q, g, g) - L 2 (q A )q 4 + C^-g 1 + g 2 + g 4 + const) (4.7) 



The Birkhoffian (|4.7|) is dissipative and not regular. 

Indeed, there exists a smooth function Eq^ : TM — > R of the form 
is, 



E o„ (<?, q) 



that 



L 1 (? 1 q 3 )(q 2 + q 3 )(dq 2 + dq 3 ) + I i 2 (g 4 )g 4 rfg 4 + / L 3 (g 2 )g 2 rfg 2 



L!(g 2 ,g 3 )rfg 2 dg 3 



C 1 (q\q 2 ,q i ){dq x -dq 2 + dq*) + J C 2 (g x ,q 3 ){dq x + dq 3 ) + 
CatfW ~ J J C[{q\ g 2 , q A ){-dq l dq 2 + dg^g 4 - dg 2 dg 4 ) - 



(4.9) 



such that (|2.16(l is satisfied with 

D = Ri{q 3 )dq 3 
Because the function R\ satisfies Q4.2[l . we obtain 

£,(g,g)g J =i?r(g 3 )g 3 >0 



(4.10) 



that is, H4.9|) is indeed a dissipative vertical 1-form. 

We are in the case where the circuit considered has one loop which contains only 
capacitors. We note that for the Birkhoffian 1)4. 7J) . the first row of the matrix 



dq* 



i,j=l,2,3,4 



contains only zeros. Therefore, det 



and the 



i,j=l,2,3,4 
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Birkhoffian 14.711 is not regular. 

Using the first relation in (|4.7I) . we now define a 3-dimensional M c C M c by 

M c = {q = (q 1 , q 2 , q 3 , q 4 ) G M c / Ci (g 1 - q 2 + q 4 + const) + 

C 2 (q 1 +q 3 + const) + C 3 {q l ) = 0Ql.ll) 

By the implicit function theorem, we obtain a local coordinate system on the 
reduced configuration space M c . Taking q 1 := q 2 . q 2 
Birkhoffian has the form u) c = Y] . ^ Qjdcp , where 



<7 3 , <j 3 := <? 4 , the 



Qi{q,q,'<l) = L x {q +q) + L 3 (q 1 ) q + L x {q + q )q - 

C 1 U(q\q\?)-q 1 +? + const) 
Q2(q,q,q) = L 1 (q 1 +q)q + L x {q + q )'q + Ri(q ) + 

C 2 (/(g\o/ 2 ,g 3 ) + g 2 + const) 
Q 3 (^,t) = L 2 (f)f + C 1 (/( g - 1 ,g 2 ,g 3 )- g - 1 + g 3 + C on S i) (4.12) 

/ : U C R 3 — ► R 1 being the unique function such that /(go) = Qo, <?o 6 R-j 
and C7i(/(«) - g 1 + g 3 + const) + C 2 (/(«) + f + const) + C 3 (f{q)) = 0, Vg = 
{q 1 ,q 2 ,q 3 ) 6 C/, with J7 a neighborhood of go = (<Zo) 9oi 9o)- 

We have shown in section 3 that in this case the reduced Birkhoffian l|4.12|l 
is dissipative and regular. 

The relations (|3.45|) are satisfied for this example, thus, there exists a function 
E 0uj (q, q) such that (|3.41|) is fulfilled, with the dissipative 1-form given by 



We calculate 



det 



dq 1 



i,j'=l,2,3 



D = Ri(q)dq 2 



Li{q + q) +i 3 (? 1 ) L^q+q) 
L 1 (q +q ) L x {q +q ) 
















(4.13) 



(4.14) 



Because L\, L 2 , L 3 : R — > R\{0}, the determinant above is different from zero, 
then, the reduced Birkhoffian given by l|4.12Jl is regular. ■ 
If the nonlinear resistor is voltage controlled, that is, 



l[i] = 9ti 

we obtain, instead of l|4.7|l . the following implicit Birkhoffian 



(4.15) 



Qi(q, q, q) = Ciiq 1 - q 2 + q 4 + const) + C 2 (q 4 + q 3 + const) + C^q 1 ) 
Q 2 (q, q, q) = [L x {q 2 + q 3 ) + L 3 (q 2 )} q 2 + L x {q 2 + q 3 )<f- 

C\ (q 1 -q 2 + q 4 + const) 
Q 3 (q, q, q) = L x {q 2 + q 3 )q 2 + L x {q 2 + q 3 )q 3 + C 2 {q l + q 3 + const) + v x (4.16) 
Q 4 {q, q, q) = L 2 (q 4 )q 4 + Ci(q 4 - q 2 + q 4 + const) 



q 3 = 9ti(«i) 
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We suppose that the nonlinear voltage controlled resistor satisfies %\.\{x)x > 0, 
for any x ^ 0. Thus, as in the case of current controlled resistors, there exists a 
function E 0ui (q, q) given by l|4.8|) , such that i|2.16|l is satisfied with the dissipative 
implicit 1-form 

D(q, (?) = Vidq 3 , 

(4.17) 

9 3 = Wi(«i) 

Therefore, the Birkhoffian (|4.16|) is dissipative. The Birkhoffian l|4.16[l is also 
not regular. ■ 

Example 2) This example is based on the following oriented connected graph 




Figure 2 



We have r = 5, fc = 2, p ~ 2, n = 5, to = 4 , 6 = 9. We choose the reference 
node to be Vq and the current directions as indicated in Figure 2. We cover the 
associated graph with the loops I\, I2, h, h- Let V = (Vi, V2, V3, V4, V5) G R 5 
be the vector of node voltage values, 1 = (im, i( a ), i Q ) € R 5 x R 2 x R 2 be the 
vector of branch current values and v — (vm, v^a),v a ) 6 R 5 x R 2 x R 2 be the 
vector of branch voltage values. 

The branches in Figure 2 are labelled as follows: the first, the second and the 
third branch are the linear resistive branches, Ru ni , Ru n2 ; Riin 3 j the fourth and 
the fifth branch are the nonlinear resistive branches R4, R5, the second, the 
sixth and the seventh branch are the inductive branches Li, L2, and the last two 
branches are the capacitor branches Ci, C2. The incidence and loop matrices, 
B G 97t 95 (R) and A e 9Jl 94 (R), write as 








-1 


1 








\ 




( 1 











\ 







1 















1 








-1 












-1 












1 


-1 
























-1 












1 















-1 


1 







A = 





1 












-1 


1 
























1 















1 












-1 


1 









1 



























1 




V 











-1 


1 


/ 




V 





1 





/ 



(4.18) 



One has rank(£>) = 5, rank(A) = 4. 

Except from the resistors in the first loop which are considered linear, all devices 
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are nonlinear and are described by the relation I3.3[l . 13. ty . i|3.6[) . We suppose 
that Ri, R2, R3 > 0, are distinct constants, C\,Ci : R — > R\{0}, Li,L 2 : R — > 
R\{0}, smooth invertible functions and i?4, -R5 : R — > R smooth functions such 
that, for any i^O 

i? 4 (a;)x > 0, i?5(a;)x > 



(4.19) 



The first set of equations i|3.9[l has the form 



-i(i) +Qi = 

+I[2] +1(1) = 
I[l] - I[3] - I[5] = 
![5] + x (2) ~ Q2 : 
-I[4] + Q2 = 



(4.20) 







The relations 1|3.12[1 . (|3.14|) read as follows for this example 

T = 1, 5, I( a ) := Q (a ) 

„6 „ JT „ „8 



X 



i[r] := 

„5 .. 



![T], 



a = 1,2 

9 .. 



Qi, x 



Q 2 



(4.21) 
(4.22) 



Q[i],..., X :=Q[ 5 ],X :=Q(i), X :=Q (2 ),ar 

Using the equations from l|4.20|l . we define the 4-dimensional afhne-linear con- 
figuration space M c . In view of the notations (|4.21l) . (|4.22l) . we integrate these 
5 equations and solving them in terms of 4 variables, we obtain, for example, 



x = x 



-x 8 +const, x 3 — x 1 —x 5 +const 1 x A — x 9 +const 1 x e — x a +const, 
x 9 + const. Thus, a coordinate system on M c is given by 

9 



q 1 :=x\q 2 :=x 5 ,q 3 :=x 8 ,q 4 



x 



(4.23) 



The matrix of constants TV from (|3.15|) is given by 



( 1 








\ 


1 





-1 





1 


-1 

















1 





1 














1 








-1 





1 








1 





V 








1 / 



Therefore, in terms of the g-coordinates (I4.23|) . 



we define the Birkhoffian uj c — Qi(q, q, q^dq 1 + Q 2 {q, q, q)dq 2 + Qs{q, q, q)dq 3 
Q 4 (q,q,q)dq 4 of M c as in (E23>, with l EHSj l. that is, 



Qi(q, q, q) = (ri + R 2 + R3M 1 - R3<? 
Q2{q,q,q) = L 2 (-q 2 



R29 



L 2 (-q 2 



•4\ -4 

q )q 



Rag 1 +R39 2 + R 5 (q 2 ) 
4 (q, q, q) = -L 2 (-q 2 + q 4 )q 2 + L 2 (-q 2 + q 4 )q 4 + R 4 (q 4 ) + C 2 (q 4 ) (4.24) 



h(q, q, q) = £l(g 3 )<f - R^ 1 + R 2 q 3 + Ci(g 3 ) 



The Birkhoffian (|4.24() is dissipative and not regular. 

Indeed, there exists a smooth function En : TM — > R of the form 



that 
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IS. 



E „(q,q) = L 1 (q 3 )q 3 dq 3 + L 2 (q\ q 4 )(-q 2 + q 4 )(-dq 2 + dq 4 ) 



L' 2 {q\q 4 ){-q 2 +q 4 Wdf+ L 2 (q 2 , q 4 )dq 2 dq 4 + 



CitfW+ / C^dq 



(4.25) 



such that (|2.1(i|) is satisfied with the dissipative 1-form defined by 

D = [Al + (q 1 - <f )R2 + (q 1 - q 2 )R 3 ] dq 1 + [-(q 1 - q 2 )R 3 + 

R 5 (q 2 )] dq 2 - (q 1 - q 3 )R 2 dq 3 + R^dq 4 (4.26) 

In view of the assumptions l|4.19|l and of Ri,R2, R3 > 0, we get 

Dj{qAW =Ri(g 1 ) 2 + R 2 (g 1 -g 3 ) 2 + R 3 (9 1 -g 2 ) 2 + i?4(g 4 )g 4 + i 1 , 5(<z 2 )'z 2 >0 (4.27) 

Therefore, the vertical 1-form in l|4.26|l is dissipative. 

We are in the case where the considered circuit has one loop which contains only 
resistors. We note that for the Birkhoffian l|4.24|) . the first row of the matrix 



dq' 



i,j=l,...,4 



contains only zeros. Therefore, dct 



— and the 



i,j'=l,...,4 



Birkhoffian i|4.24[l is not regular. 
Using the first relation in (|4.24|) . we now define M c C M c by 

M c = {q= (q^q 2 , q 3 ,q 4 ) G M c j (r x + R 2 + R-^q 1 - R 3 <Z 2 - R 2 <7 3 + Cj = 0} (4.28) 

where ci is a real constant. 

On the reduced configuration space M c , in the coordinate system given by q 1 :— 
q 2 , q 2 := q 3 , q 3 :— q 4 , the Birkhoffian has the form cD c = Qidq 1 + Q 2 dq 2 + Q 3 dq 3 

Ql{q,'q,l) = L 2 (-q + q)q ~ L 2 (-q + q)q ' + (<t 2 + €3) q - <t 3 q ' + R 5 (q) 
Q 2 {qAA) = L x (q 2 )f -<L 3 £ + {£l + U)q +Ci{q 2 ) 



• 1 -.3, -;3 



3 (q, q,q) = -L 2 (-q + q )q + L 2 (-q + q )q + R 4 (q ) + C 2 (q 3 ) 



(4.29) 



R1R2 



R1R3 



R1+R2+R3 ' 



R1+R2+R3 



where we denote the constants £1 
R 2 R3 

Ri+R 2 +R 3 ' n-im 

As we have stated in section 3, the Birkhoffian given by l|4.29|l is still dissipa- 
tive. The function Eq^ (q, q) has the same form (I4.25f) written in the coordinates 
q. The relation (|3.t)2l) is satisfied with the dissipative 1-form defined by 



D 



q£ 2 + (q - q)£a + dq 1 + \q<t x - (q - q)£ 3 ] dq 2 + 



Ri(q )df 



(4.30) 



2N 



The vertical 1-form above is dissipative, as can be seen as follows: For Ri > 0, 
R 2 > 0, R 3 > 0, we get £i > 0, € 2 > 0, C 3 > and together with (|4~T3|) yield 



Dj(q,q)q = ^(q f + € 3 (q - 'q ? + ^(q ) 2 + R 5 (q )q + R 4 (q )q >0 (4.31) 



The Birkhofhan given by (|4.29() is not regular, since the determinant 



dct 



dQj 
dq 



z,j=l,2,3 



L 2 (q) ^0 

Li(q) 
-L 2 (q) 



-L 2 {q) 
^ 
L 2 (q) 



(4.32) 



This result does not come unexpected because the considered network has also 
a loop which contains only resistors and capacitors, formed by R4, C 2 , R5, Ru n3 . 
In order to regularize the Birkhofiian l|4.29|l . we introduce an inductor in series 
into this loop, described by the following relation between the current and the 
voltage: v = £1(1)^, C\ : R — * R\{0} being smooth invertible function. 
This means that this loop will have one more node and one more branch. The 
number of branches for the graph associated to the circuit, increases by one, 
that is, there will be b = 10 branches, and the number of nodes increases by one 
as well, that is, we will have n — 6. But the cardinality of a selection of loops 
which cover the whole graph remains m = 4. After the calculation we arrive at 
the reduced configuration space defined by (|4.28|l . On the reduced configuration 
space M c , in the coordinate system given by q 1 :— q 2 , q 2 :— q 3 , q 3 := q 4 , the 
Birkhofiian ui^. xt has the components Qi(q, q, q), Q 2 (q, q, q) given by (|4.29J) and 
the expression of Qs(q, q, q) becomes 



Q3(q,q, q) 



We now calculate 



det 



dq 



-L 2 (-q + <f )^ + \Ci(q) + L 2 (-q + q) 



R^q °) + C 2 (q 3 ) 



L 2 (q) ^ ~L 2 (q) 
--'■■< -L 2 {q) Cx(q)+L 2 (q) 



q 



(4.33) 



(4.34) 



Because Ci, L\, L 2 : R — > R\{0} the determinant above is different from zero, 
then, the Birkhofiian Cj^. xt is regular. ■ 
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